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Abstract. We present a model for nonlocal diffusion with Neumann boundary condi- 
tions in a bounded smooth domain prescribing the flux through the boundary. We study 
the limit of this family of nonlocal diffusion operators when a rescaling parameter related 
to the kernel of the nonlocal operator goes to zero. We prove that the solutions of this 
family of problems converge to a solution of the heat equation with Neumann boundary 
conditions. 



1. Introduction 

The purpose of this article is to show that the solutions of the usual Neumann boundary 
value problem for the heat equation can be approximated by solutions of a sequence of 
nonlocal "Neumann" boundary value problems. 

Let J : M. N -» R be a nonnegative, radial, continuous function with J RN J(z)dz = 1. 
Assume also that J is strictly positive in -8(0, d) and vanishes in WL N \ -8(0, d). Nonlocal 
evolution equations of the form 

(1.1) u t (x,t) = (J * u — u)(x,t) = / J(x — y)u(y,t) dy — u(x,t), 

and variations of it, have been recently widely used to model diffusion processes. More 
precisely, as stated in [TU], if u(x,t) is thought of as a density at the point x at time t 
and J(x — y) is thought of as the probability distribution of jumping from location y to 
location x, then J RN J(y — x)u(y, t) dy = (J * u)(x, t) is the rate at which individuals are 
arriving at position x from all other places and — u(x, t) = — J RN J(y — x)u(x, t) dy is the 
rate at which they are leaving location x to travel to all other sites. This consideration, in 
the absence of external or internal sources, leads immediately to the fact that the density 
u satisfies equation (jl.lj) . For recent references on nonlocal diffusion see, [T], [2], 0, jl], 
0, 0, 0, PH, [12], PH, [H and references therein. 

Given a bounded, connected and smooth domain Q, one of the most common boundary 
conditions that has been imposed in the literature to the heat equation, u t = Am, is the 
Neumann boundary condition, du/drj(x,t) = g(x,t), x e dQ, which leads to the following 
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classical problem, 



(1.2) 



u t — Am = 
du 




in Qx(0,T), 
on dQx(0,T) 



u(x, 0) = u (x) 



in fi. 



In this article we propose a nonlocal "Neumann" boundary value problem, namely 



where G(x, £) is smooth and compactly supported in £ uniformly in x. 

In this model the first integral takes into account the diffusion inside Q. In fact, as 
we have explained, the integral f J(x — y)(u(y,t) — u(x,t))dy takes into account the 
individuals arriving or leaving position x from or to other places. Since we are integrating 
in Q, we are imposing that diffusion takes place only in Q. The last term takes into 
account the prescribed flux of individuals that enter or leave the domain. 

The nonlocal Neumann model (jl.Hj) and the Neumann problem for the heat equation 
(jl.2|) share many properties. For example, a comparison principle holds for both equations 
when G is nonnegative and the asymptotic behavior of their solutions as t — > oo is similar, 
see 0. 

Existence and uniqueness of solutions of (jl.Hj) with general G is proved by a fixed point 
argument in Section 2. Also, a comparison principle when G > is proved in that section. 

Our main goal is to show that the Neumann problem for the heat equation (|1.2jl can 
be approximated by suitable nonlocal Neumann problems (jl.Hj) . 

More precisely, for given J and G we consider the rescaled kernels 



z JB(0,d) 

which is a normalizing constant in order to obtain the Laplacian in the limit instead of a 
multiple of it. Then, we consider the solution u £ (x,t) to 




(1.4) 




with 






(1.5) 



< 



u e (x,0) = u {x). 



We prove in this paper that 

u e — > u 

in different topologies according to two different choices of the kernel G. 
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Let us give an heuristic idea in one space dimension, with Q = (0, 1), of why the scaling 
involved in (jl.4|) is the correct one. We assume that 

J(y)ydy 



poo pu pi 

/ G(l,l-y)dy=- / G(0,-y)dy= / 

Jl J-oo JO 



and, as stated above, G(x, •) has compact support independent of x. In this case (jl.5j) 
reads 

If 1 If 

u t (x,t)= — J e (x-y)(u(y,t)-u(x,t))dy + - G £ (x,x - y) g(y,t) dy 



1 



+oo 



+ - / G £ (x, x-y) g(y, t) dy := A £ u(x, t). 
8 Ji 

If x G (0, 1) a Taylor expansion gives that for any fixed smooth u and e small enough, 
the right hand side A £ u in ()1.5|) becomes 

l r 1 

A £ u(x) = — J £ (x-y) (u(y) - u(x)) dy « u xx (x) 
£ Jo 

and if x = and e small, 

A £ u(0) = - 2 [ J £ (-y) (u{y) - u(0)) dy + - [ G £ (0, -y) g{y) dy w ^(n x .(0) - p(0)). 

Analogously, A £ u(l) ~ (C2/e)(— u x (l) + fl'(l))- However, the proofs of our results are 
much more involved than simple Taylor expansions due to the fact that for each e > 
there are points x G f2 for which the ball in which integration takes place, B(x,de), is 
not contained in Q. Moreover, when working in several space dimensions, one has to take 
into account the geometry of the domain. 

Our first result deals with homogeneous boundary conditions, this is, g = 0. 

Theorem 1.1. Assume g = 0. Let Q be a bounded C 2+a domain for some < a < 1. 
Let u G C 2+aA+a/2 (U x [0,T]) be the solution to flL2j) and let u £ be the solution to (fPf 
wii/i J £ as above. Then, 

sup ||w e (-,£) - w(-,£)||l->(q) -> 
te[o,T] 

as e — » 0. 

Note that this result holds for every G since (7 = 0, and that the assumed regularity in 
u is standard if wo G C 2+a (fi) and duo/dr] = 0. See, for instance, [IT] . 

We will prove Theorem II. II bv constructing adequate super and subsolutions and then 
using comparison arguments to get bounds for the difference u £ — u. 

Now we will make explicit the functions G we will deal with in order to consider g ^ 0. 

To define the first one let us introduce some notation. As before, let Q be a bounded 
C 2+a domain. For x G Vt £ := {x G Q \ dist(x, dtt) < de} and e small enough we write 
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x = x — s dr]{x) where x is the orthogonal projection of x on dQ, < s < e and r](x) is 
the unit exterior normal to Q at x. Under these assumptions we define 

(1.6) G 1 {x,£) = -J{£)r ] (x)-{i forxeQ £ . 

Notice that the last integral in (jl.5j) only involves points x G Q £ since when y ^ Q, 
x — y G supp J e implies that x G £l e . Hence the above definition makes sense for e small. 

For this choice of the kernel, G = G\, we have the following result. 

Theorem 1.2. Let Q be a bounded C 2+a domain, g G C 1+a ^ 1+a)/2 (dtt x [0,T]) ; u G 
C 2+a ' 1+a/2 (tt x [0,T]) the solution to JOJ), /or some < a < 1. Let J as before and 
G(x,£) = Gi(x,£), where G\ is defined by (ll.fjjl . Le^ u e 6e i/ie solution to (ll.5jl . Then, 

sup ||M e (-,t) - u(-,t)||z,i(n) -> 

ie[0,T] 

as £ — >• 0. 

Observe that Gi may fail to be nonnegative and hence a comparison principle may not 
hold. However, in this case our proof of convergence to the solution of the heat equation 
does not rely on comparison arguments for (|1.3|) . If we want a nonnegative kernel G, in 
order to have a comparison principle, we can modify (Gi) E by taking instead 

(G a ) 6 (a;, £) = (Gi) e (x, + «eJ e (0 = \ MO [-v{x) ■ £ + ke 2 ) . 

Note that for x G Q and y eR N \Q, (Gi) £ (x, x - y) = ~M X ~ V) ( _? ?(^) " i x ~ V) + K£2 ) 
is nonnegative for e small if we choose the constant k as a bound for the curvature of dQ, 
since \x — y\ < de. As will be seen in Remark |4.1[ Theorem II . 21 remains valid with (G\) £ 
replaced by (Gi) £ . 

Finally, the other "Neumann" kernel we propose is 

G(x,0 = G 2 (x,0 = C 2 J(0i 

where C 2 is such that 

(1.7) / / J(z) (C 2 - z N ) dz ds = 0. 

JO J{z N >s} 

This choice of G is natural since we are considering a flux with a jumping probability that 
is a scalar multiple of the same jumping probability that moves things in the interior of 
the domain, J. 

Several properties of solutions to ()1.3|) have been recently investigated in [S] in the case 
G = G 2 for different choices of g. 

For the case of G 2 we can still prove convergence but in a weaker sense. 

Theorem 1.3. Let Q be a bounded C 2+a domain, g G C 1+a ^ 1+a)/2 (dtt x [0,T]) ; u G 
C 2+a ' 1+a/2 (Tt x [0,T]) the solution to (fT^jl . for some < a < 1. Let J as before and 
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G(x, £) = G 2 (x, £) = C 2 J(£), where C 2 is defined by ()1.7|) . Let u e be the solution to (JT3J). 
T/ien, /or eac/i i € [0, T] 

u e (x, t) — ^ it (a;, t) * —weakly in L°°(fl) 

as e — > 0. 

The rest of the paper is organized as follows: in Section|2]we prove existence, uniqueness 
and a comparison principle for our nonlocal equation. In SectionElwe prove the uniform 
convergence when g = 0. In Section 0] we deal with the case G = G\ and finally in 
SectionEJwe prove our result when G = G 2 . 



2. Existence and uniqueness 

In this section we deal with existence and uniqueness of solutions of (jl.3j) . Our result 
is valid in a general L 1 setting. 

Theorem 2.1. Let Q be a bounded domain. Let J e L 1 (M N ) and G E x R N ). 

For every u G L 1 ^) and g G L^ c ([0, oo); L 1 (M 7V \ fl)) there exists a unique solution u of 
(jl.3l) such that u G C([0, oo); L l {VL)) and u(x, 0) = Uq(x). 

As in j2j and |Hj , existence and uniqueness will be a consequence of Banach's fixed point 
theorem. We follow closely the ideas of those works in our proof, so we will only outline 
the main arguments. Fix to > and consider the Banach space 

x tQ = c([o,t );L 1 M) 

with the norm 



\w\ 



o<t<t |L1 ^ 



We will obtain the solution as a fixed point of the operator T U0i9 : X to —* X to defined by 
T UQ>g (w)(x,t) = u (x) + [ [ J(x-y)(w(y,s)-w(x,s))dyds 

G(x,x -y)g(y,t) dyds. 



o Jn 

t 



(2.1) 

+ 

h Jr n \q 

The following lemma is the main ingredient in the proof of existence. 

Lemma 2.1. Let J and G as in Theorem, PHI Let g, h G L°°((0,t o ); L 1 (M N \ ft)) and 
Uq, Vq G L 1 (J1). There exists a constant C depending only on Q, J and G such that for 
w,z G X to , 

(2.2) \\\T UOig (w)-T VOth (z)\\\ < \\u -v \\ L i +Ct (\\\w - z\\\ + \\g - h\\ L ^ i{0M)) . L i {R N\ n)) ) . 
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Proof. We have 



\T UOjg (w)(x,t) - T VOth (z)(x,t)\ dx < / \u (x) - v (x)\ dx 



+ 



+ 




o Jn 

t 



J{x-y) (w(y,s) - z(y,s)) - (w(x,s) - z(x,s)) dyds 



dx 



\G(x,x -y)\\g(y,s) - h(y, s)\ dydsdx. 



'n Jo 

Therefore, we obtain (|2.2j) . 



□ 



Proof of Theorem \2.1\ Let T = T uo g . We check first that T maps X tQ into X to . From 
fl2~T]) we see that for < t x < t 2 < t , 



\T(tc)(U) - :/ ? f»')f-/- J .)|| i i (n) <A I 2 / \w(y,s)\dyds + B [* 

h Jn Jt! 



\g(y,s)\dyds. 



On the other hand, again from (j2.1J) 

\\T(w)(t) - U \\ L i (n) < Ct{\\\w\\\ + \\g\\L°°((0,t );Li(RK\n))}- 

These two estimates give that T{w) G C([0,t ]; L l (Vt)). Hence T maps X tQ into X to . 

Choose to such that Cto < 1. From Lemma 12.11 we get that T is a strict contraction 
in X to and the existence and uniqueness part of the theorem follows from Banach's fixed 
point theorem in the interval [0, to]- To extend the solution to [0, oo) we may take as 
initial datum u(x,t ) G and obtain a solution in [0,2t ]- Iterating this procedure 

we get a solution defined in [0, oo). □ 

Our next aim is to prove a comparison principle for (jl.3j) when J, G > 0. To this end 
we define what we understand by sub and supersolutions. 

Definition 2.1. A function u G C([0, T); is a supersolution of (jl.3|) ifu(x,0) > 

u (x) and 

u t (x,t)> / J(x -y)(u(y,t) -u(x,t)) dy+ / G(x,x - y)g(y,t)dy. 



Subsolutions are defined analogously by reversing the inequalities. 

Lemma 2.2. Let J, G > 0, uq > and g > 0. If u G C(Q x [0, T]) is a supersolution to 
fOjl . then u>0. 



Proof. Assume that u(x,t) is negative somewhere. Let v(x,t) = u(x,t) + et with e so 
small such that v is still negative somewhere. Then, if we take (xo,to) a point where v 
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attains its negative minimum, there holds that to > and 

v t {x ,t ) = u t (x ,t ) + e > / J(x -y)(u(y,t ) -u(x ,t ))dy 

Ja 

J(x- y)(v(y,t ) - v(x ,t )) dy > 
which is a contradiction. Thus, u > 0. □ 



Corollary 2.1. Let J, G > and bounded. Let uq and Vq in L X (Q) with uq > no and 
g, h E L oo ((0,T);L 1 (R N \ Q)) with g > h. Let u be a solution of ((EHJ) wift im'fca/ 
condition Uq and flux g and v be a solution of (ll.Hjl nnt/i initial condition vq and flux h. 
Then, 

u > v a.e. 

Proof. Let w = u — v. Then, w is a supersolution with initial datum no — Vq > and 
boundary datum g — h > 0. Using the continuity of solutions with respect to the initial 
and Neumann data (Lemma EH} and the fact that J G L°°(R N ), G G x R N ) we 

may assume that u, v G C(Q x [0, T]). By Lemma f2. 21 we obtain that to = it — v > 0. So 
the corollary is proved. □ 

Corollary 2.2. Let J, G > and bounded. Let u G C(fi x [0,T]) (Vesp. v) be a 
supersolution (resp. subsolution) of 1)1. Hj) . Then, u > v. 

Proof. It follows the lines of the proof of the previous corollary. □ 

3. Uniform convergence in the case y = 

In order to prove Theorem II .11 we set w £ = u £ — u and let u be a C 2 + a ' 1 + a / 2 extension 
of u to M w x [0, T]. We define 



L £ (n) = ^ / J E (x-y)(v(y,t) -v(x,t))dy 



and 



L £ (n) = I Je{x - y)(v{y,t) - v(x,t))dy. 



Then 



tu t £ = L £ (u e ) - An + - / G e (x, x - y)g(y, t) dy 



e _ 

L £ {w e ) + L £ {u) - Au + i / G e (x, a; - n)y(n, t) dy 
1 



e 2 



- y) (u(y, t) - u(x, t)) dy. 



Or 



w £ t -L £ {w £ ) =F £ {x,t), 
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where, noting that Am = Am in Q, 

F e (x, t) = L £ (u) - Au + - / G £ (x, x - y)g(y, t) dy 

e Ju N \n 

Je( x ~ y) (u(y, t) - u(x, t)) dy. 



1 



R N \n 



Our main task in order to prove the uniform convergence result is to get bounds on F £ . 

First, we observe that it is well known that by the choice of C\, the fact that J is 
radially symmetric and u G C 2+a,1+a / 2 (R JV x [0,T]), we have that 

(3.1) sup \\L £ (u) - Au\\ Lao(n) = 0(e a ) . 

te[o,T] 



In fact, 



' ^ j( ) (u(y,t) — u(x,t)) dy — Au(x,t) 



becomes, under the change variables z = (x — y)/e, 
Ci 



/ J (z) (u(x — ez, t) — u(x, t)) dy — Au(x, t) 

and hence (jH.ljl follows by a simple Taylor expansion. 

Next, we will estimate the last integral in F e . We remark that the next lemma is valid 
for any smooth function, not only for a solution to the heat equation. 

99 

Lemma 3.1. If 6 is a C 2+a ' 1+a ^ 2 function on R N x \0,T] and — = h on d£l, then for 

ar\ 

x G £l £ = {z G Q | dist (z, dQ) < de} and e small, 

- [ J e {x - y) (6{y, t) - 9(x, t)) dy = - [ J £ {x - y) V {x) ■ ^t—^-h{x, t) dy 

e JR N \n £ 

D&9 f _ ^\({y — ( (x - x)^p- 

1/91=2 

where x is the orthogonal projection of x on the boundary of Q so that \\x — y\\ < 2de. 

Proof. Since 9 G C 2+a > 1+a / 2 (R N x [0,7]) we have 

9(y, t) - 9{x, t) = 9(y, t) - 9(x, t) - (9(x, t) - 9{x, t)) 

D@9 

= W9(x, t)-{y-x)+J2 -2~& *) [(V - *Y ~(x- xf] 

|/3|=2 

2+a\ i n/IU „,I|2+Q\ 



f t / \ F) f3 9 , . r, (y — x) x a / (x — x) x b! , , 
+ / J e {x-y) V — ~{x,t) -) -(" ") dy + 0{e c 



+ 0(\\x - x\\ 2+a ) + 0(\\x-y 
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s 2 



1 f J e {x-y) (0(y, t) - 9(x, t)) dy = - [ J £ (x - y)V9(x, t) ■ dy 
Jw N \n £ Jr n \q. £ 

+ / J £ (x-y)J2 C-^—^T-C-—^) dy + 0(e a ). 

JM. N \n fzr„ lie s J 



'R N \n 

d6 

Fix x G Q e . Let us take a new coordinate system such that nix) = ex- Since — = h 

Of] 

on <9Q, we get 



/ J e (x-y)V9{x,t).^—El d y 

JR N \n £ 



= [ J £ (x-y)r ] (x). { y-^l h{ x,t)dy+ [ J £ (x - y) V 9 Xi (x, t)^—^- dy. 
Jr n \q, £ JR N \n ~[ 

We will estimate this last integral. Since f2 is a C 2+a domain we can chose vectors e 1: 
e 2 , 6 TV— 1 so that there exists k > and constants such that 



N-l 



£ 2< fe(z) n \y N - (x N + ^ - ^) 2 ) < -/«r 2+Q j c a 



i=i 



Therefore 

N-l 



[ J £ (x-y)(f^e Xt (x,t) { y^-)d !J 
jR N \n K ~l e / 

r N_1 ( - ) 

/ Je{x -y)(J2 e *i( 2 > *) ) d y 

V(RJv\n)n||i/ JV -(* JV +i;^ 1 /i(x)(vi-*i) 2 )|<«e a+a } V i= i 5 7 

/■ JV_1 ( _ ^ 



= /l + / 2 . 

If we take z = (y — x)/e as a new variable, recalling that xjy — x N = es, we obtain 

JV-l . 

J I Ujv-(«+ei;£i 1 /i(S)(«i) 2 )|<«e 1+a j 



i=l 
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On the other hand, 

N-l „ 

h = CiY] 6xi (x, t) J (z) Zi dz. 

i= i J[z N -{s+eY.t\ 1 /.(i)( 2 ,) 2 )>^+°) 

Fix 1 < % < N — 1. Then, since J is radially symmetric, J(z) Zi is an odd function of 
the variable z^ and, since the set |zjv — (s + £ Yl!i=i fi(.%)(. z i) 2 ) > K£ 1+Q? j is symmetric 
in that variable we get 

h = 0. 

Collecting the previous estimates the lemma is proved. □ 

We will also need the following inequality. 
Lemma 3.2. There exist K > and e > such that, for e < e, 

(3.2) / J £ (x - y)rj(x) ■ dy > K [ J £ (x - y) dy. 

Jr n \q £ Jr n \q 

Proof. Let us put the origin at the point x and take a coordinate system such that 
7)(x) = e N . Then, x = (0, — /i) with < ji < de. Then, arguing as before, 



/ J £ (x - y)rj(x) ■ — — —dy=[ J £ (x-y) V '' 
Jr n \q, £ Jr n \q, 



dy 



/ J £ (x - y) VN + - dy + / J £ (x-y)- 

J{y N >Ke 2 } £ JR N \nn{\y N \<KE 2 } 



e 



dy 



> j J £ (x-y) y -^dy-Ce. 

'{y N >Ke 2 } £ 



Fix Ci small such that 

1 /' 

2 



/ J(z) z N dz > 2ci / J(z) dz. 

JUn>0\ il0<^«<2ci> 



'{z N >0} J {0<z N <2c 1 } 

We divide our arguments into two cases according to whether \i < c\e or /i > c\S. 
Case I Assume jj, < c\E. In this case we have, 

/ J £ (x - y) VN + - dy = d / J(z)z N dz 

J{y N >Ke 2 } £ J{z n >ks+!±} 

(3.3) = C\ / J{z)z N dz— I J{z)z N dz 

\J{z N >0} J{0<Z N <KS+^} / 

> Ci ( j J{z)z N dz — 2c\ j J(z) dz J > f J{z)z^dz. 

\J{z N >0} J{0<z N <2 Cl } J 2 J{ ZN>0 } 
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Then, 



if e is small enough and 



(v — x) 

J £ (x - y)r](x) dy-K I J E {x-y) dy 



J(z) z N dz- K) -Ce>0 



{z N >0} 



K < 



4 



J(z) z N dz. 



{z N >0} 



Case II Assume that fi > C\E. For y in \QC\ B(x, de) we have 

yN 



Then, 



> — KE. 



Je( x ~ y) VN + ~ dy-K I J £ (x - y) dy 



r n \u 

> (ci — ke) 



J £ (x — y)dy — K / J £ (x - y) dy 
R N \n 



(d-KE-K) j J £ (x-y)dy>0, 
JR N \n 



if e is small and 



K < 



Cl 



This ends the proof of (|3.2|) . 
We now prove Theorem 11.11 



□ 



Proof of Theorem \l.l\ We will use a comparison argument. First, let us look for a super- 
solution. Let us pick an auxiliary function v as a solution to 



'v t - Av = h(x,t) in fix(fl,T), 
dv 

^v(x, 0) = Vx(x) 



on dttx(0,T), 
in Q. 



for some smooth functions h(x,t) > 1, gi(x,t) > 1 and Vi(x) > such that the resulting 
v has an extension v belongs to C 2+a ' 1+a ^ 2 (M N x [0,T]), and let M be an upper bound 
for v in Q x [0,T]. Then, 



«t = L £ t> + (Av — L £ v) + 



J £ {x - y) (v(y, t) - v(x, t)) dy + h(x, t). 
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Since Av = Av in Q, we have that v is a solution to 

v t - L £ v = H(x,t,e) in ftx(0,T), 
v(x, 0) = f i(x) in $7, 

where by Lemma 13. II and the fact that fa > 1, 

1 



H(x,t,e) =(Av - L £ v) 



Je(x - y) (v(y, t) - v(x, t))dy + h(x, t) 



J £ (x - y)r]{x) ■ — — —gi(x, t) dy 



|/3|=2 



J £ {x-y)}^—{x,t) (— ^) -(— 



)1 ^) 



> 



0lOM) 



J £ (x - y)rj(x) ■ dy - D 1 



J £ (x - y) dy 



l-Ce c 
1 



for some constant D\ if e is small so that Ce a < 1/2. 

Now, observe that Lemma 13.21 implies that for every constant Cq > there exists Eq 
such that, 



(y - x) 



J £ (x-y)r)(x) ■ 

R N \n £ 



dy-Co / J £ (x-y)dy > 0, 

R N \n 



if e < £ . 

Now, since g = 0, by (|3.1J) and Lemma f3. II we obtain 



\F £ \ < Ce a + 



J £ (x -y)}^ -z-{x,t) {——-) - {-^—) 



dy 



101=2 



jR N \n 



Given 8 > 0, let v$ = Sv. Then v$ verifies 

(v s ) t - L e v s = 5H(x,t,e) in Qx(0,T), 
v§(x, 0) = Svi(x) in Q. 

By our previous estimates, there exists e = e (5) such that for e < e , 

\F £ \ < 5H(x,t,e). 
So, by the comparison principle for any e < Eq it holds that 

-M<5 < -va <w £ <v 5 < M5. 

Therefore, for every 5 > 0, 

— M5 < liminf w £ < limsupw e < MS. 

" e^O 
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and the theorem is proved. □ 

4. Convergence in L 1 in the case 67 = G\ 

First we prove that F £ goes to zero as e goes to zero. 
Lemma 4.1. If G = G\ then 

F £ (x,t)^0 in L^dO^^L 1 ^)) 

as e — > 0. 

Proof. As G = G x = - J(0 7/(5) ■ f, for a; e fi 6 , by (jSHJ) and Lemma EH 

F £ (x, t) =- / J e {x- y)v(x) (g(y, t) - g(x, t)) dy 

£ JR N \n e 



- Je{x-y) > — -(x,t) ( ) -( ) dy + 0{e ). 

Jr n \q |^T 2 i y e 6 J 

As g is smooth, we have that F £ is bounded in Q £ . Recalling the fact that |fl e | = 0(e) 
and F £ (x,t) = 0(e a ) on Q \ Q £ we get the convergence result. □ 

We are now ready to prove Theorem 11.21 

Proof of Theorem M.HX In the case G = G\ we have proven in Lemma [4. II that F £ — >• in 
L x (n x [0, T]). On the other hand, we have that w e = u e — u is a solution to 

w t - L £ (w) = F £ 

w(x,0) = 0. 



Let z e be a solution to 



Then — z s is a solution to 



By comparison we have that 



z t -L £ (z) = \F £ \ 
z(x, 0) = 0. 

z t -L £ (z) = -\F £ 
z(x,0) = 0. 



-z £ <w £ <z £ and z £ > 0. 
Integrating the equation for z £ we get 

\\z £ (-,t)\\ L i(Q) = / z £ (x,t)dx= / / \F £ (x,s) \ dsdx . 
Jn Jn Jo 



sup \\z £ (-,t)\\ L i {n) -> 
te[o,r] 



Applying Lemma 14.11 we get 
as £ — > 0. So the theorem is proved. □ 
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Remark 4.1. Notice that if we consider a kernel G which is a modification of G\ of the 
form 

G £ (x,0 = (G 1 ) £ (x,Z)+A(x,£,e) 

with 

\A(x, x — y, e) \ dy — > 

in as e — > 0, then the conclusion of Theorem is still valid. In particular, we 

can take A(x, £, e) = keJ £ (C,). 

5. Weak convergence in L 1 in the case G = G 2 
First, we prove that in this case F £ goes to zero as measures. 
Lemma 5.1. If G = G 2 then there exists a constant C independent of e such that 

T 

F e (x, s) \ dxds < C. 



o Jn 
Moreover, 

F £ (x,t) ^ as measures 
as e — > 0. That is, for any continuous function 9, it holds that 

[ [ F £ (x,t)9(x,t)dxdt -> 
Jo Jn 

as e — > 0. 



Proof. As G = G 2 = C 2 J(£) and g and u are smooth, taking again the coordinate system 
of Lemma 13.11 we obtain 

1 f t t \(n / ±\ Vn - xn 



F £ (x, t) =- / J £ (x- y) (c 2 g(y, t) - y*—E» g fa t)) 



-f J E (x -y)^u x .(x,t)— — —dy 
£ JR N \n ~[ e 



ux - y )j2 uk^ly _ ( fc^M dy + 0(£ c 



f J £ (x-y)(c 2 g(x,t)- — — —g(x,t) 
R N \n v £ 

AT-1 



£ 



- f J £ ( x -y)Tu x Xx,t) { -^^-dy + 0(l) X n £ + 0(£ a ). 
£ Ju N \n ~[ £ 



NONLOCAL DIFFUSION PROBLEMS 15 

Let 

B E (x, t) := / J £ (x - y) [C 2 g(x, t) - — — —g(x, t) ) 



/ J £ (x~y) y2u Xi (x,t)— — — dy. 
Jn*\n ~l £ 



Proceeding in a similar way as in the proof of Lemma f3. II we get for e small, 

Je( x - y) (c 2 g(x, t) - Vn - Xn g(x, t)j 



9(x,t) [ Ux-y) (c 2 - {m - XN) ) dy 

J {m N \n)n{\y N -x N \<KE 2 } \ £ J 

+ 9frt) [ J £ {x-y) (c 2 - iyN - XN) ) dy 

J (R N \Q)n{y N -x N >0} \ 



£ 



uu-.t) I •;-(•*• -/;)K' 2 - (y * Xn) ) <hi 

' (W. N \Cl)n{Q<y N -x N <K,e 2 } 



And 



C ig (x,t) / J(z)(C 2 - z N )dz + 0(e) X n E 

'{ZN>S} 



[ J £ (x-y) S^u Xi (x,t)— — —dy 

n-i . / 

= ^2u Xi (x,t) / J e (x-y)- 

i=l J {\vn-xn\<ke 2 } 



e 



N-l . 

+ ^u Xl (x,t) / Je(x-y) 

J {y N -x N >Ks 2 } 



{Vi ~ Xi) 



i=i 

N-l 



e 



dy 



= Ci ^2u Xi (x,t) / J(z)zidz + 0(e)xn e 

j = l J {Z N -S>K£} 

= I 2 + 0(e) X n £ . 
As in Lemma f3. II we have I 2 = 0. Therefore, 

B £ (x,t) = dg(x,t) / J(z)(C 2 -z N )dz + 0(e) X n e . 

J{z N >s} 

Now, we observe that B £ is bounded and supported in Q £ . Hence 

/ [ \F £ {x,r)\dxdT <- [ [ \B £ (x,r)\dxdT + Ct\Q £ \+Ct\Q\e a <C. 
Jo Jn £ Jo Jn E 

This proves the first assertion of the lemma. 
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Now, let us write for a point i£Sl £ 

x = x — /ir)(x) with < /i < de. 

For e small and < fi < de, let dS^ be the area element of {x G f2 1 dist (x, dQ) 
Then, dS^ = dS + 0(e), where dS is the area element of dQ. 



So that, taking now /i = se we get for any continuous test function 6, 
B £ (x,t)9(x,t)dxdt 



1 <- T 



o Jn, 



0(e) + d [ [ g(x,t)9(x,t) [ [ J(z)(C 2 - z N ) dz ds dS dt 
Jo Jan Jo J{z N >s} 

= 0(e) ^ as e -> 0, 



since we have chosen C 2 so that 



/ / J(z) (C 2 - z N ) dz ds = 0. 

JO J{z N >s} 



Now, with all these estimates, we go back to F £ . We have 

F e (x,t) = l -B £ (x,t) + 0(l) X n £ + 0(e a ). 

Thus, we obtain 

/ / F £ (x,t)6(x,t)dxdt -> ase^O. 
Jo Jn £ 

Now, if <j(r) is the modulus of continuity of 0, 



F £ (x,t)6(x,t)dxdt = [ [ F £ (x,t)6(x,t)dxdt 

Jo Jq e 

+ [ [ F £ (x,t){6(x,t) -6(x,t))dxdt 
Jo Jq £ 

< f [ F £ (x,t)9(x,t)dxdt + Ca(e) [ [ \F £ (x,t)\ dx dt -> as e -> 
Jo Jn e Jo Jn e 

Finally, the observation that F £ = 0(e a ) in £7 \ f2 e gives 

/ / F £ (x,t)9(x,t)dxdt -> ase^O 
Jo Jn\n e 

and this ends the proof. 

Now we prove that w e is uniformly bounded when G = G 2 . 
Lemma 5.2. Let G = G 2 . There exists a constant C independent of e such that 



p £ |lz,°°(f}x[0,T]) ^ C- 
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Proof. Again we will use a comparison argument. Let us look for a supersolution. Pick 
an auxiliary function v as a solution to 

'v t ~ Aw = h(x,t) in fix(0,T), 



(5.1) 



dv 



gi{x,t) 



., on dOx(O.T), 

or] 

^v(x, 0) = Vi(x) in fi. 

for some smooth functions h(x,t) > 1, t>i(a;) > Mo(x) and 

2 _ 

9i{x,t) > — (C 2 + l^max |^(a;,*)| + 1 (If as in (Q ) 

such that the resulting t> has an extension 5 that belongs to c ,2+a ' 1+Q! / 2 (R JV x [0,T]) and 
let M be an upper bound for v in Q x [0, T]. As before u is a solution to 

v t - L £ v = H(x,t,e) in ftx(0,T), 
f (x, 0) = f i(x) in fi, 



where if verifies 



H(x,t,e) > 



gi(x,t) 



(y — gr\ 

J e (x - y)f]{x) dy - D 1 



R N \n 



J £ (x-y)dy) + \. 



So that, by Lemma \'A.2\ 



H(x,t,e) > ( 



g 1 (x,t)K 



J £ (x -y)dy+- 



for e < e. 

Let us recall that 



F e (x,t) = L £ (u) - Au + 



9l 

e 

1 



Je{% - y)g(y,t) dy 



Je( x ~ v) (u(y, t) - u(x, t)) dy. 



Then, proceeding once again as in Lemma f3. II we have, 



\F £ (x,t)\ < 



\g{x,t)\C 2 



< 



J £ (x -y)dy + 

+ Ce a + C [ J £ {x-y)dy 
(C 2 + 1) 



9(x,t)\ 



Je{x - y)\n{x) 



{y-x) 



dy 



S 9Hx[0,T] 

if e < e, by our choice of g\. 



max \g(x,t)\+C / J £ (x — y) dy + Ce c 



R N \n 



J £ {x-y) dy + Ce a 
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Therefore, for every e small enough, we obtain 

\F e (x,t)\ <H(x,t,e), 
and, by a comparison argument, we conclude that 

—M < -v(x,t) < u £ (x,t) < v(x,t) < M, 
for every (x,t) G Q x [0, T\. This ends the proof. 

Finally, we prove our last result, Theorem 11.31 

Proof of Theorem M.iA By Lemma f5. II we have that 

F £ (x,t) ^ as measures in fix [0, T] 

as e —* 0. 

Assume first that ^ G Cg +Q (fi) and let y5 e be the solution to 

w t — L £ w = 
w(x, 0) = 4>(x). 



□ 



Let be a solution to 



' (p t — Acp — 
< ^ = 

or] 

k y(x,0) = 



Then, by Theorem ll.il we know that ip £ —>■ (p uniformly in f2 x [0,T]. 
For a fixed t > set <p E (x, s) = (p £ (x, t — s). Then yj e satisfies 

ip s + L e (/) = 0, for s < t, 

<p(x,t) = i){x). 



Analogously, set (p(x, s) 



<p(x,t — s). Then cp satisfies 
f (p t + A<p = 

or) 

K (p(x,t) = ip(x). 
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Then, for w £ = u £ — u we have 

/"* f dw £ f l f d^ £ 

w £ (x, t) ip(x) dx = / / — — (x, s) <£ £ (x, s) dx ds + / / -jr—(x, s) w £ (x, s) dxds 
<» Jo Jn V s Jo Jn V s 

L £ (w £ )(x, s)<f e (x, s) dx ds + / / F e (x,s) ip e (x,s) dxds 



o Jn 

t 



o Jn 



<9y? £ 
ds 



[x, s) w e (x, s) dx ds 



o Jn 

L E ((p e )(x, s)w £ (x, s) dxds + / / F e (x, s) (p £ (x, s) dxds 



o Jn 

t 



o Jn 



+ 



dps 



. -(x, s) w e (x, s) dx ds 
o Jn ds 

F e (x, s)ip e (x, s) dx ds. 



o Jn 



Now we observe that, by the Lemma f5. 11 



F E (x, s)(p £ (x, s) dx ds 



o Jn 



< 



F e (x, s)ip(x, s) dx ds 



o Jn 



+ sup \\<p e (x,s)-(p(x,a)\\ L «>(ti) 

o<s<t jo jq 



F £ (x, s) \ dxds — > 



as e — > 0. This proves the result when ip G C^+^n;). 

Now we deal with the general case. Let ip G Choose ip n G Cq +01 (H) such that 

ip n — > ip in We have 



w e (x, t) ijj(x) dx 



< 



W £ (x, t) 1pn{x) dx 



\L 1 {n)\\W \\L°°{0)- 



By Lemma f5. 21 {w £ } is uniformly bounded, and hence the result follows. 



□ 
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